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ABSTRACT 

Neutron star (NS) merger events can eject considerable amounts of neutron-rich matter, whose 
expansion allows heavy nuclei up to uranium and thorium to be formed by the rapid neutron-capture 
process. The radioactive decay of the products of this r-process heats the ejecta material and allows it 
to become potentially observable as a source of thermal electromagnetic radiation. Here we investigate 
for the first time systematically the dynamical mass ejection during and shortly after the NS collision 
in dependence on the uncertain properties of the nuclear equation of state (EoS) by employing a set of 
40 representative, microphysical high-density EoSs in relativistic, hydrodynamical merger simulations. 
The NS compactness, which is a characteristic EoS property and can be expressed by the radius -R1.35 
of nonrotating NSs of 1.35 M Q , turns out to be a crucial parameter that determines the ejecta mass. 
NSs with smaller radii R1.35 ("soft" EoS) collide more violently and eject systematically higher masses. 
These range from ^10~ 3 M Q to ~10 -2 M Q for symmetric 1.35-1.35 Mq binaries with -R1.35 between 
16 km and 11km, and from ^5 x 10 -3 Af© to -2 x 1O -2 M for asymmetric 1.2-1.5 Af© binaries. 
Correspondingly, the bolometric peak luminosities of the optical transients of symmetric (asymmetric) 
mergers vary between 3 x 10 41 ergs -1 and 14 x 10 41 ergs -1 (9 x 10 41 ergs -1 and 14.5 x 10 41 ergs -1 ) and 
are reached with effective spectral temperatures of 13. 000-19. 000K on typical timescales between ^2 h 
and ~12h. If these electromagnetic signals with absolute bolometric magnitudes from —15.0 to —16.7 
can be measured, the tight correlation of their properties with those of the merging NSs might provide 
valuable constraints on the high-density EoS. The wide range of ejecta velocities and kinetic energies 
in our models (between 15% and 50% of the speed of light and between ^5 x 10 49 erg and 10 51 erg, 
respectively) suggests variations by factors of several 100 for the peak of the radio flux expected from 
the ejecta interaction with the circumstellar medium. In contrast, despite differences in the ejecta 
dynamics and thermodynamics, the r-process nucleosynthesis exhibits a remarkable robustness. In all 
investigated cases more than 93% of the ejecta mass are converted to neutron-rich nuclei with nuclear 
mass numbers A >130 and a uniform, nearly solar abundance pattern. By the r-process content of 
the Galaxy and the average production per event the Galactic merger rate is limited to 4 x 10 -5 yr -1 
(4 x 10 -4 yr -1 ) for a soft (stiff) NS EoS, if NS mergers are the main source of heavy r-nuclei. The 
production ratio of radioactive 232 Th to 238 U, which we determine for the first time in the discussed 
context, attains a stable value of 1.64-1.67, which does not exclude NS mergers as potential sources 
of heavy r-material in the most metal-poor stars. 

Subject headings: equation of state — hydrodynamics — nuclear reactions, nucleosynthesis, abun- 
dances — stars: abundances — stars: neutron 



1. INTRODUCTION 

Neutron star (NS) merger events are among 
*_h . the most promising candidates for the first di- 
& • rect measurement of a gravitational-wave signal 
with the upcoming Advanced LIGO and V IRGO 
interferometric instruments (lAcernese et all 120061 : 
IHarrv fe LIGO Scientific Collaborationll2010|) . and they 
are considered as likely origin of short gamma-ray 
bursts and their afterglows as a conseque nce of ultra- 
relativistic, collim ated outflows (see e.g. iNakarl [20071 : 
iKann et "all 1201 ID . Moreover, they are possible sources 
of different kinds of electromagnetic signals in the 
precursor of the merging and in its aftermath as a 
consequence of magnetohydrodynamical effects, magne- 
tospheric interactio ns, relativistic matter outflows , or NS 
crust phenomena dLipunov fc Panchenkol 119961 : IVietril 
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produced by hot ejecta gas, for 
potent i ally o bservable opti- 
transients (ILi fc Paczvhskil 119981: [Kulkarnil 120051: 
IMetzger et "all 120101: IMetzger fc Bergerl l2012l Tand the 
interaction of the ejecta cloud with the circumstellar 
medium is expected to create ra dio flares that might 
be detectable for periods of year s (INakar fc Pirad l20lH: 
iPiran et alj|2012t iRosswog et alj|2012f ). Observations of 
such signals could help to pinpoint the exact celestial 
locations of NS mergers (thus, e.g., supporting the 
analysis of data taken by gravitational- wave detectors), 
and repeated measurements of signals that can be 
unambiguously linked to NS mergers would help to 
constrain the still highly uncertain rate of such events 
in the local universe. 
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During the merging of two NSs a small fraction of 
the system mass, typically 0.1-1 per cent, can be- 
come gravitationally unbound and can be ejected on 
the dynamical timescale of mil l iseconds (IRu ffcrt et al. 
19971: IRosswog et alj 11991 [2000t IRuffert fe Jankal 120011: 
Oechslin et al.l I2007t IRosswog et al.l I2012t IPiran et al.l 



2JH2J; [Rosswog |2012t EptokezakleLalJ |2013j). Because 



such material is likely to possess a high neutron ex- 
cess, it has been proposed as a possible site for the 
creation of the heaviest, neutron-rich elements, which 
are form ed by the rapid neutron capture proc ess (r- 
process) (jLattimer et al.lll977t lEichler et al.lll989f ) (sim- 
ilarly, also NS-black hole mergers were suggested as 
sourc es of r-process matter (jLattimer fc Schrammlll974l 
119761) 1 The radioactive decay of these freshly syn- 
thesized r-process nuclei should heat the ejecta and 
thus lead to an optical transient (|Li fc Paczvhski| 11998 



also needs to be understood which fraction of the ejecta 
is robustly converted to r-process material and whether 
the final abundances are always compatible with the solar 
element distribution, which agrees amazingly well with 
the r-process abundance pattern in metal-poor stars fo r 
atomic numbers Z ~ 55-90 (see e.g. lSneden et al.ll2008| ). 

Newtonian as well as relativistic studies showed that 
the mass ratio has a significant effect on the amount 
of matter that can become unbound (IJanka et all 
1991 IRosswog et al.lll99l 120001: IRuffert fc Jankal 12001 



Oechslin et al 
20111 



2007: 



[Roberts et al 
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IPiran et all 120121: IRosswog et al.l 1201 
Korobkin et al.l l2012t IRosswog! 120121: iHotokezaka et al.l 
20131) . Such investigations, however, were performed 



only with a few exemplary models fo r high-density 
matt e r in NSs (IRosswog et al.1 120001: lOechslin et ail 
[20071 IGorielv et al.l 120111: IHotokezaka et al.l I2013D or 
Kulkarnil 1200a IMetzger et al.H2010l: IRoberts et al.ll20Hl: even only with a single NS EoS (IRoberts et all l20lTt 
Goriely et all 1201 ID . The properties of such events de- Piran et all 120121: IRosswog et al.l I2012t iKorobkhT et al.l 



pend on the fraction of the material that can be con- 
verted to radioactive species. Moreover, the peak lumi- 
nosity, the timescale to reach the emission peak, and the 
effective temperature at the radiation maximum, as well 
as the radio brightness that accompanies the deceleration 
of the expelled gas during its coasting in the stellar envi- 
ronment, depend sensitively on the ejecta mass and ex- 
pansion velocity. Detailed hydrodynamical merger mod- 
els are needed to calculate these quantities and to de- 
termine the nucleosynthesis conditions in the unbound 
material. 

Concerning their role as sources of heavy elements bi- 
nary NSs collisions have recently moved into the focus of 
interest because the astrophysical sources of the r-process 
elements have not been identified yet and core-collapse 
supernova simulations continue to be unable to yield the 
extreme conditions for forming the heaviest neutron- 
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However, despite this promising situation a variety 
of aspects need to be clarified before the question can 
be answered whether NS mergers are a major source 
or even the dominant source of heavy r-process ele- 
ments. On the one hand the merger rate and its evo- 
lution during the Galactic history are still sub j ect to 
considerable uncertainties (see e.g.. lAbadie et al.l (|2010f ) 
for a compilation of recent estimates), and it is unclear 
whether NS mergers can explain the early enrichment of 
the Galaxy by r-process elements as observed in metal- 
deficient stars (jArgast et al.l 120041 ). On the other hand 
it remains to be determined how much mass is ejected 
in merger events depending on the binary parameters 
and, in particular, depending on the incompletely known 
properties of the equation of state (EoS) of NS matter. It 



nuclear EoS for a quantitative assessment of the dy- 
namical mass ejec tion can be conclud ed from published 
calculations (e.g. IGorielv et al.l 120111 ). These calcula- 
tions, however, also suggest that the nuclear abundance 
pattern produced by r-processing in the ejecta may be 
largely insensitive to variations of the conditions in the 
ejecta. 

It is important to note that quantitatively reliable 
information on the ejecta masses and their dependence 
on the binary and EoS properties require general 
relativistic (GR) s imulations. Newt o nian results 

in the literature (IRosswog et all [l99l IJanka et all 

1991 IRuffert fc Jankal 120011: IRoberts et al.l 12011 



Korobkin et all l2012t IRosswog et all 120121 : IPiran et all 
20121 : IRosswog! 120121 ) exhibit significant quantitative 
and qual itative differen c es co m pared to rela t ivistic 
models (lOechslin et all 120071: IGorielv et all I20ITI: 
IHotokezaka et al.ll2013l) . Newtonian calc ulations tend to 
overesti mate the ejec ta masses in general (IRosswog et all 
19991: IJanka etal.1 [19991: IRuffert fc Jankal 120011: 
Roberts et alJl201lUKorobkin et al 1120121 IRosswog et all 
20121 : IPiran et all 120121 : IRosswog! 120121 ). This can be 



understood because of several facts. First, the structure 
of NSs in GR is considerably more compact than 
that of Newtonian stars. For instance, a NS with a 
gravit ational mass of 1.35 M w described by the LS220 
EoS (jLattimer fc Swestvl 119911 ) possesses a circumfer- 
ential radius of 12.6 km, whereas the corresponding 
Newtonian star has 14.5 km. Second, GR gravity is 
stronger and the merging of two NSs is therefore more 
violent. The difference can be expressed in terms of the 
gravitational binding energy of a nuclcon on the surface 
of the considered 1.35 M Q NSs, which is ~200MeV 
in the GR case compared to only <~130 MeV for the 
Newtonian model. Third, GR forces merger remnants 
beyond a mass limit to collapse to black holes on a 
dynamical timescale. Such an effect cannot be tracked 
by Newtonian models. These differences arc of direct 
relevance for the collision dynamics and the possibility 
to unbind matter from the inner and outer crust regions 
of the merging NSs. 

It is the purpose of this paper to explore the in- 
fluence of the high-density EoS on the ejecta proper- 
ties in a systematic way, i.e., we will determine ejecta 
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masses and the nucleosynthesis outcome for a large set 
of NS matter models, applying them in relativistic NS 
merger simulations. Most of th ese EoSs were already em- 
ployed in our previou s works (jBauswein fc Jankalf20Tl 
iBauswein et al.1120121 ). They were chosen such that they 
provide as completely as possible a coverage of the pos- 
sibilities for NS properties (expressed by corresponding 
mass-radius-relations) which are compatible with present 
ob servational const r aints (e.g. , the 1.97 M© NS discovery 
of iDemorest et al.l (120101)) and theoretical understand- 
ing (lLattimer k Prakashll2010L [20071 iSteiner et al.ll201Ct 
iHebeler et al.ll2010D . In our study we will focus on sym- 
metric 1.35-1.35 M© systems and will compare them with 
asymmetric 1.2-1 .5 mergers. Beca use population 
synthesis m ode ls (Belc zvnski et a l. 2008) a nd pulsar ob- 
serva tions (jThorsett k Chakrabartvl 119991 : iZhang et al.1 
120 111 ) suggest that the double NS population is strongly 
dominated by systems of nearly equal-mass stars of about 
1.35 Mq each, the average NS merger event can be well 
represented by a 1.35-1.35 M configuration, and a clari- 
fication of the EoS dependence of ejecta masses, r-process 
yields, and properties of electromagnetic counterparts of 
NS mergers seems to be more important than a wide vari- 
ation of binary parameters. Nevertheless, we will also 
present results of a more extended survey of binary mass 
ratios and total masses for some representative EoSs. 

In our work we will exclusively concentrate on NS-NS 
mergers, but the discussed phenomena sh ould play a 
role a lso for NS-black hole coalescence dJanka et al.1 
19991 iLeel 



Faber et al 



Piran et al 



20001: IRosswog et al.l 



Foucart et al.l 



NS 



mergers 



(lEast et al 



Rosswog et al.l I2012D and eccentric 



- 12012 MEast k Preform! [20121 
IRosswog et al.l l2012t IRosswog 12012ft . However, while 
the existence of double NS systems is established 
by observations, progenitors of NS-black hole and 
eccentric NS mergers have not been observed yet 
and the rates of such types of events are even more 
uncertain than those of coalescing binary NSs. In 
investigating the latter we will only consider the phase 
of dynamical mass ejection between about the time 
when the two NSs collide until a few milliseconds 
later. During this phase hydrodynamical and tidal 
forces (shock compression, pressure forces, gravitational 
interaction) are responsible for the mass shedding of 
the merging objects. Once the remnant has formed, 
however, differential rotation i s expected to strongly 
amplify the magnetic fields (e.g. iPrice k Rosswogil2006l: 
Anderson et al . 2008; L iu et alll2008t iGiacomazzo et al.l 



20TI " and viscous energy dissipation is likely to provide 



additional heating, enhancing the neutrino emission that 
accompanies t he secular evolut i on of t he post-merger 
confi g uration (IRuffert k Jankal 119991: iSetiawan et all 
l200i IDessart et al.l 120091: IShibata k Sekiguchil I2012D . 
As a consequence the merger remnant will experience 
mass loss due to neutrino energy deposi t ion in the near- 
surfa ce regions (IRuffert k Jankal Il999t [Setiawa.n et all 



2001 IDessart et al.l 120091 : iWanajo fc Jankal 120121: 



Shibata fe Sekiguchil 120121 ) (similar to the neutrino 



driven wind of proto-neutron stars emerging from stellar 
core collapse) and due to magnctohydrodynamical 
outflows. Both mechanisms will add ejecta to the mass 
stripped during the dynamical interaction of the system 
components, but the details of the secular evolution 



and the associated mass loss will be very sensitive to 
the EoS-dcpendcnt stability properties of the merger 
remnant, i.e., to the qu estion whether th e remn ant is 
a hypermassive NS (see iBaumgarte et all ((2000) for a 
definition) or whether and when it collapses to a black 
hole-torus system. These questions lie beyond the scope 
of the present work. 

Our paper is organized as follows. In Sect. [2] a brief 
summary of the numerical methods and microphysics 
ingredients of our NS merger simulations is given. In 
Sect. [3] we present our results for the relation between 
dynamical mass loss and NS (EoS) properties, provide 
a detailed description of the mass-loss dynamics in our 
relativistic models (drawing comparisons to Newtonian 
results), discuss the influence of an approximate treat- 
ment of thermal effects in the EoS, and evaluate the 
mass ejection for three selected, representative EoSs in 
merger simulations for a wider space of binary masses 
and mass ratios in order to determine the population- 
integrated mass loss. In Sect. [4] we describe results of 
nuclear network calculations performed for a subset of 
our merger models and draw conclusions on the Galac- 
tic merger rate and the production of long-lived radioac- 
tive species ( 232 Th, 235 U, 238 U) used for stellar nuclc- 
ocosmochronomctry. Finally, we present values for the 
heating efficiency of the merger ejecta by radioactive de- 
cays of the nucleosynthesis products and apply them in 
Sect.[5]to estimate the properties (peak luminosity, peak 
timcscalc, effective temperature at the maximum lumi- 
nosity) of the optical transients that can be expected 
from the expanding merger debris. We also briefly dis- 
cuss the implications of our simulations for radio flares. 
Finally, a summary and conclusions follow in Sect [6] 

2. NUMERICAL MODEL AND EQUATIONS OF STATE 

The simulations of our study are performed with 
a relativistic Smooth Particle Hydrodynamics (SPH) 
code, i.e. the hydrodyna mical equations are evolved 
in a Lagrangian man ner (jOechslin et al.l [20021 I2007T : 
IBauswein et al.ll2010al ). The Einstein field equations are 
solv ed imposing conformal f latness of the spat ial met- 
ric (jlsenberg fe Nested Il980l: iWilson et all 11996ft . and a 
gravitational-wave backreaction scheme is used to ac- 
count for energy and angular mo mentum losses by th e 
emission of gravitational radiation (jOechslin et aLll2007t ) . 
The code evolves the conserved rest-mass density p* . the 
conserved specific momentum Ui, and the conserved en- 
ergy density r, whose definitions evolve the metric poten- 
tials and the "primitive" hydrodynamical quantities, i.e. 
the rest-mass density p, the coordinate velocity Vi, and 
the specific internal energy e. The system of relativis- 
tic hydrodynamical equations is closed by an EoS, which 
relates the pressure P = P(p,T,Y c ) and the specific in- 
ternal energy e = e(p, T, Y e ) to the rest-mass density p, 
the temperature T and the electron fraction Y c . The 
temperature is obtained by inverting the specific inter- 
nal energy e = e(p, T, Y e ) for given p and Y c . Changes of 
the electron fraction ar e assumed to be slow compared to 
the dynamics (see e.g. IRuffert et al.lll997l) . and the ini- 
tial electron fraction, which is defined by the neutrinoless 
beta-equilibrium of cold NSs, is advected according to 
^[f = (j| defines the Lagrangian, i.e. comoving, time 
derivative) . 

The EoS of NS matter is only incompletely known 
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and numerical studies rely on theoretical prescriptions 
of high-density matter. This work surveys a repre- 
sentative sample of 40 microphysical EoSs, which have 
been derived within different theoretical frameworks 
and make different assumptions about the composition 
of high-density matter and the description of nuclear 
int eractions. Mo s t of t he employed EoSs are listed 
in iBauswein et al.l ([2012D . where details can be found, 
while some new models are introduced below. Because 
of the one-to-one correspondence between the EoS and 
the mass-radius relation of nonrotating NSs, it is conve- 
nient to characterize EoSs by the resulting stellar prop- 
erties. Stellar quantities as integral properties of an 
EoS are in particular useful to classify the dynamics 
of NS mergers and the accompanying gravitational- wav e 
signals (jBauswein fc Ja nka 2012L IBauswein et alll2012T ). 
For this reason we will adopt the same approach also 
for this investigation. Considering for instance NSs with 
a gravitational mass of 1.35 Mq, the stellar radii R1.35 
vary from 10.13 km to 15.74 km for the different EoSs 
of our sample. The maximum mass M max of nonrotat- 
ing NSs obtained for these EoSs ranges from 1.79 Mq 
to 3.00 Mq. In terms of their stellar properties the em- 
ployed EoSs of our study show a large variation (seethe 
mass-radius relations in Fig. 4 of IBauswein et all (|2012f) V 
Note that we do not apply any selection procedure for 
choosing the EoSs, except that we require a maximum 
mass above w 1.8 Mq. This limit is chosen because of 
the firm discovery of a pulsar with a gravi tational mass 
of (1.97 ± 0.04) Mq (jDemorest et al.|[20Tot) . EoSs which 
yield a maximum mass below this limit are practically ex- 
cluded by this observation. Nevertheless, we accept them 
(at least down to M max sa 1.8 Mq) for our investigation 
because we expect that at densities relevant in a typical 
NS merger these models stil l provide a viable desc ription 
of high-density matter (see [Bauswein et al. 2 0121) . Note 
that c ompared to our previous study in IBauswein et al.l 
poll we extend our EoS survey by including also 
the models TM1, TMA, NL3, DD2, SFHO and SFHX 
of IHempel fe Schaffner-Bielichl (|2010D IHempel et all 



(|2012D and ISteiner et al.l (120121), rely i ng on the inter- 
action s described in Sugahara fc Tqkil (119941) iToki et all 
JT995I). lLalazissis et all (1199711 ifybef et al.l (|2010D 



and ISteiner et al.r ~( 2012D . Moreover, w e inclu de the 
BSk20 and BSk21 EoSs of IGorielv et all (1201(1 . The 
maximum masses M max resulting for these EoSs are 
2.21 Mq, 2.02 Mq, 2.79 Mq, 2.42 Mq, 2.06 Mq, 
2.13 M , 2.16 M Q and 2.28 Mq (order as listed above), 
while the radii of cold 1.35 Mq NSs are 14.49 km, 
13.86 km, 14.75 km, 13.21 km, 11.92 km, 11.98 km, 
11.74 km and 12.54 km, respe c tively. From our sam- 
ple of EoSs in IBauswein et all (|2012l) we do not con- 
sider the SKA EoS (because of its restriction to densities 
above 1.7 x 10 9 g/cm 3 ) and EoSs which are not compati - 
ble with the pulsar observation of IDemorest et al.l (|2010f ) 
and directly form a black hole after merging. We also ex- 
cl ude absolutely stable s trange quark matter. We refer 
to IBauswein et all ([2009D for the particular implications 
of ejecta from strange quark star mergers. 

Only 12 out of the considered 40 EoSs describe ther- 
mal effects consistently and provide the dependence of 
thermodynamical quantities on the temperature and the 
electron fraction. Instead, the majority of models consid- 



ers matter at zero temperature and in equilibrium with 
respect to weak interactions (i.e. for beta-equilibrium for 
neutrino- less conditions). Because temperature effects 
become important during the merging of the binary com- 
ponents and during the subsequent evolution, we employ 
an approximate treatment of thermal effects for those 
EoSs which are given as barotropic relations. This pro- 
cedure supplements the pressure by an additional ideal- 
gas component to mimic thermal pressure support, and 
it requires to choose a corresponding ideal-gas index r t h- 
Appropriate values for r t ,h are in the range o f 1.5 to 2 
for high-density matter (Bauswei n et al.lfeOlOal ). The un- 
certainties connected to the use of this approximate tem- 
perature descript ion and the choice of the ideal-gas index 
were examined in lBauswein et all ()2010al ). where also de- 
tails about the exact implementation c an be found. 

Fro m population synthesis studies (jBelczvnski et all 
120081 ) and in agreement w it h pu ls ar observa- 
tions (jThorsett fc Chakrabartvl Il999t iZhang et all 
120111) it is expected that binaries with two NSs with 
gravitational masses of about M\ w M2 1.35 Mq are 
the most abundant systems. For this reason we focus 
in our EoS survey on such equal-mass binaries, albeit 
we also explore the influence of a system asymmetry by 
considering 1.2-1.5 Mq binaries. For a selected subset of 
EoS models the full range of possible binary parameters 
is investigated, varying the single component masses 
from 1.2 Mq to approximately the maximum mass of 
NSs. 

Because of energy and angular momentum losses by 
gravitational radiation the orbits of NS binaries shrink 
and the binary components merge after an inspiral pe- 
riod, which last s roughly 100 to 1000 Myrs for the 
known systems (jLorimerl 120081 ) . The typical outcome 
of the coalescence of a 1.35-1.35 Mq binary system is 
the formation of a differentially rotating object, poten- 
tially a hypermassive NS (i.e. a NS that is more mas- 
sive than the maximum-mass rigid-rotation configura- 
tion and that is stabilized temp orarily by differential 
rotation (jBaumgarte et all l2000f )) . The merger rem- 
nant is surrounded by an extended halo structure of 
low-density material. Only four EoSs of our sample 
lead to the prompt formation of a black hole within 
about one millisecond after the collision because the rem- 
nant can not be supported against the gravitational col- 
lapse. For a description of the general dynamics and 
a more thorough discus s ion o f the collapse beh a vior we 
refe r to lOechslin et all (1200711 IBauswein et~aTI (|2010al ) 
and IBauswein et al.l (|2012l ). In this paper only initially 
nonrotating NSs are investigated because viscosity is too 
low to yield tidally locked systems. The stars in NS bi- 
naries are therefore expected to rotate slowly in compar- 
ison to the orbital angular velo city, justifying the use of 
an irrotational velocity profile (jBildsten fc Cutlerl 119921 : 
IKochaneldll99l . 

In this study we analyze the material which becomes 
gravitationally unbound during or right after merging. 
In order to estimate whether a given fluid element, i.e. 
an SPH particle, can escape to infinity, we consider 



1 



^stationary 



1 > 



(i) 



with the coordinate velocity v % , the conserved momcn- 
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turn &i, and the time-component of the eigen- velocity 
u° (in geometrical units). This expression can be de- 
rived from the hydrodynamical equations by neglect- 
ing pressure forces an d assuming a stationary met- 
ric (jOechslin et al.ll2002l ). The quantity e s tationary is con- 
served ( deata ^° nar y = and at infinity it reduces to the 

Newtonian expression for the total energy of a fluid ele- 
ment. Hence, a particle with e s tationary > will be un- 
bound. Equation ([T]) is evaluated in a time-dependent 
way and SPH particles that fulfill this criterion 10 ms 
after merging are considered as ultimately gravitation- 
ally unbound. Note that our simulations neglect a pos- 
sible (smaller) contribution to the ejecta by neutrino- 
driven winds or magnetically driven outf lows from the 
secular evolution of the mer ger remnant (|Dessart et al.l 
[20091: IWanaio fc ,Tankal[20Tl . 

3. EJECTA MASSES 

In the following we employ Eq. (fTJ to determine the 
unbound material in different merger simulations. After 
a first, steep rise of the ejecta mass shortly after the merg- 
ing of the two NSs, the mass fulfilling the ejecta criterion 
remains approximately constant (Fig. |T|). A few models, 
however, show a continuing, slow increase of the ejecta 
mass also at later times. The ejecta masses discussed 
below are computed 10 ms after merging. 

In order to determine the influence of the high-density 
EoS on the ejecta we employ approximately the same nu- 
merical resolution of about 350,000 SPH particles for all 
simulations. By using nonuniform SPH particle masses 
to model the stellar profile (more massive particles in 
the high-density core and lighter particles in the outer 
low-density layers) it is possible to achieve a better res- 
olution of low-density regions. This results in an effec- 
tive mass resolution of about 2 x 10~ 6 M , which is 
comparable to a SPH simulation of about one million 
equal-mass particles. The influence of the numerical res- 
olution is investigated by performing additional simula- 
tions with higher SPH particle numbers. For the TM1 
EoS the ejecta masses are found to be (in solar masses) 
1.67 x 10~ 3 ,1.80 x 10~ 3 ,1.71 x 10~ 3 ,2.43 x 10~ 3 , and 
2.07 x 10~ 3 for calculations with about 339 x 10 3 , 550 x 
10 3 , 782 x 10 3 , 1007 x 10 3 , and 1272 x 10 3 SPH particles. 
Determining the unbound matter 6 ms after merging for 
the APR EoS, we find ejecta masses (in solar masses) 
of 5.93 x 10" 3 , 6.08 x 10~ 3 , 6.54 x 10~ 3 , 6.69 x 10~ 3 , and 
6.14 x 10~ 3 in simulations with 339 x 10 3 , 592 x 10 3 , 782 x 
10 3 , 1007 x 10 3 , and 1272 x 10 3 SPH particles. Thus, the 
numerical resolution has an effect on the level of some 
10 per cent. This, however, is smaller than the impact 
of the EoS (see below), which is the focus of this paper. 
The nonmonotonic variations of the ejecta mass with in- 
creasing resolution indicate that statistical fluctuations 
have some influence on the ejected particle population 
as well. 

3.1. Origin of the ejecta and comparison with other 
calculations 

As can be seen in Fig. [5] most of the ejecta originate 
from the contact interface between the colliding binary 
components, which get deformed into drop-like shapes 
prior to the merging. For the 1.35-1.35 Mq binary the 
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Fig. 1. — Evolution of the minimum lapse function a (dashed 
line) and the amount of unbound matter (solid line) for the sym- 
metric 1.35-1.35 Mq merger with the soft SFHO EoS (top panel), 
the intermediate DD2 EoS (middle panel), and the stiff NL3 EoS 
(bottom panel). 
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ejecta in the shear interface between the stars are sepa- 
rate into two components, each being fed (nearly) sym- 
metrically by material from both colliding stars (top 
right panel and middle left panel). The matter in the 
cusps of the stars essentially keeps its direction of mo- 
tion towards the companion, whereas the backward part 
of the contact interface mixes with some of the compan- 
ion matter (top right panel). Both lumps of ejecta are 
squeezed out from the contact interface and expand on 
the retral side of the respective companion star, par- 
tially slipping over it (middle panels). The bulk mat- 
ter of the binary components forms a rotating double- 
core structure where the two dense cores oscillate against 
each other (not visible because o f the logarithmic density 
scal e; see e.g. the description s in lBauswein et al.l ()2010bl ) 
and iStergioulas et all (|2011[ ) ) . A first bunch of the mat- 
ter which was squeezed out from the contact sheet gets 
unbound in a first expansion phase of the rotating dou- 
ble core structure, which pushes the ejecta outward. This 
can be seen in Fig. [2] (middle right panel and bottom left 
panel) and also in the evolution of the minimum lapse 
function a (Fig. [lj , which is a measure for the compact- 
ness of the central object. As the cores separate from 
each other and the lapse function grows out of its mini- 
mum, the ejecta mass increases. A second expansion of 
the double cores unbinds a smaller amount of matter. Fi- 
nally, about two milliseconds after the first contact, the 
triaxial deformation grows into two spiral-arm like ex- 
tensions reaching out from the central remnant. These 
expand into the surrounding, low-density halo fed from 
the contact interface, push it away and unbind additional 
matter (bottom right panel of Fig.[2]and second mass-loss 
episode visible in Fig. [J). For different EoSs the different 
dynamical mechanisms contribute to the ejecta produc- 
tion with different relative strengths. For soft EoSs (top 
panel of Fig. Q} the first steep rise of the ejecta mass due 
to the expanding double core is much more pronounced, 
whereas for stiff EoSs (bottom panel of Fig. [TJ the first 
increase of the ejecta is very moderate and the late spi- 
ral arms unbind most of the ejecta. This can be seen in 
Fig.Q]for the SFHO EoS representing a soft EoS, for the 
NL3 EoS as a stiff example, and the intermediate case of 
the DD2 EoS. 

A significantly smaller fraction of the ejecta (typically 
below 25 per cent) stems from the outer faces of the merg- 
ing stars opposite to the contact layer (SPH particles at 
the outer left and outer right ends of the stellar body 
in the top right panel of Fig. [3]). During merging this 
matter at the rear of the star (SPH particles with nearly 
horizontal velocity vectors at the top and bottom merger 
tails in the middle left panel) lags behind the rotation of 
the star's center and is hit and ablated by the "nose" of 
the companion shortly after the snapshot shown in the 
middle left panel (see velocity arrows of particles with 
opposite color at the tips of the noses). This material 
gets mostly unbound in the first expansion phase of the 
oscillating double-core structure. Such type of ejecta is 
less abundant for stiff EoSs. 

Different from relativistic simulations, Newtonian 
models find the ejec ta originating mostly from the tips 
of tidal tails (see e.g. lKorobkin et aT1l2012f ) . in particular 
also in the case of symmetric binaries. Relativistic cal- 
culat i ons (within the CF C framework) ( Occhsli n et al.l 
l2007t iGorielv et al.l 120111) yield the dominant ejection 



from the contact inter face as descri b ed ab ove (see also 
the inset of Fig. 1 in IGorielv et ail (120111) V Recently 
the fu lly relativistic simulations of IHotokezaka et alj 
(|2013[) have provided further support for the ejecta ori- 
gin fro m the contact interfa ce, c onfirming the conclu- 
sions of IQechslin eTaT] (|2007h and IGorielv etaLl (pOPl . 
This points to qualitative differences between the Newto- 
nian and relativistic mass-loss dynamics with the impor- 
tant difference that in relativistic simulations all of the 
ejecta are shock-heated while in Newtonian calculations 
the cold, tidally stripped material dominates. 

A quantitative c o mpari so n between Newto- 
nian (IKorobkin et al.1 120121 IRosswog et al.1 I2012t 
iPiran et al l l2012t IR osswoe) 120 121) and relativist i c sim- 
ulations ~(|Oechslin et al.l 120071: IGorielv et all l20Tfl) 
also reveals considerable discrepancies. For instance, 
simulations of a 1.4-1.4 M Q merger with the Shen EoS 
in Newt onian theory p r oduce more than 10~ 2 Mq 
eiecta (IKorobkin et all 120121: IRosswog etHI I2012t 
IPiran et al.ll2012t IR osswo gf l20l"l . whereas the relativis- 
tic c alculations of this stu dy and in lOechslin et al.l (|2007f ) 
and IGorielv et all (|2011l) yield only a few 10" 3 M of 
unbound material for the 1.35-1.35 Mq binary with the 
same EoS. Comparing the results o f our study with 
the lik ewise relativistic calculations in IHot okezak a et al.1 
(120131) shows very good agr eement for all four EoSs used 
in IHotokezaka et al.l (|2013l ). For example, for the APR 
EoS with r t h = 2 both groups find about 5 x 10 -3 M 
of unbound matter. This is remarkable because the im- 
plementations differ with respect to the hydrodynamics 
scheme, which is an SPH (smooth particle hydrodynam- 
ics) algorithm here but a grid-based , high -resolution 
central scheme in IHotokezaka et all ((2013) . (Note 
that we employ the confo rmal flatness app r oxima tion 
whereas the calculations in IHotokezaka et al.l (|2013|) are 
conducted within full general relativity.) These findings 
provide confidence in the results on the quantitative 
level and point towards fundamental differences between 
Newtonian and relativistic treatments. Such differences 
are not unexpected because NSs are more compact 
in general relativity than in Newtonian gravity. The 
stronger gravitational attraction prevents the formation 
of pronounced tidal tails at the outer faces of the 
colliding stars and increases the strength of the collision. 



3.2. Equation of state dependence 

Several NS EoSs have been employed in merger simu- 
lations by different groups, but a large, systematic inves- 
tigation of the EoS dependence of the ejecta production 
is still missing in particular with a consistent description 
of thermal effects. For a given EoS the radius -R1.35 of 
a nonrotating NS with 1.35 Mq is a characteristic quan- 
tity specifying the compactness of NSs. Therefore, we 
use -R1.35 to describe the influence of the high-density 
EoS on the amount of NS merger ejecta. 

The upper left panel of Fig. [3] displays the amount of 
unbound material as a function of R1.35 for all 40 EoSs 
used in our study. Red crosses identify EoSs which pro- 
vide the full temperature dependence. The black sym- 
bols correspond to barotropic zero-temperature EoSs, 
which arc supplemented by a thermal ideal-gas compo- 
nent choosing T t h = 2 (sec Sect. [2]). Results based on 
the same zero-temperature EoS but with T t h = 1-5 are 
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Fig. 2. — Merger and mass ejection dynamics of the 1.35-1.35 Mq binary with the DD2 EoS, visualized by the color-coded conserved 
rest-mass density (logarithmically plotted in g/cm 3 ) in the equatorial plane. The dots mark SPH particles which represent ultimately 
gravitationally unbound matter. Their positions are projections of the three-dimensional locations anywhere in the merging stars onto the 
orbital plane. Black and white indicate the origin from the one or the other NS. For every tenth particle the coordinate velocity is indicated 
by an arrow with a length proportional to the absolute value of the velocity (the speed of light corresponds to a line length of 50 km). The 
time is indicated below the color bar of each panel. Note that the side length of the bottom panels is enlarged. 
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given in blue at the same radius -R1.35. Small symbols 
indicate results for EoSs which are excluded by the pul- 
sar mass measurement of lDemorest et ail (|2010f) . Circles 
mark cases which lead to the prompt collapse to a black 
hole. 

One can recognize a clear EoS dependence of the ejecta 
mass, where EoSs with a high compactness of the NSs 
lead to an enhanced production of unbound material. 
The ejecta mass can be as big as about 0.01 Mq for 
symmetric mergers with a total binary mass M to t = 
Mi + M 2 = 2.7 M Q . EoSs with relatively large NS 
radii lead to outflow masses of about 0.001 to 0.002 Mq. 
For EoSs with approximately the same -R1.35 the ejecta 
masses show a scatter of up to 0.003 Mq. However, 
considering only EoSs with a fully consistent descrip- 
tion of thermal effects (red symbols) the variations are 
smaller. Only one simplified EoS (eosAU) leads to a 
prompt collapse of the merger remnant and yields signif- 
icantly smaller ejecta masses (circles). Using the radius 
Ri.e of a nonrotating NS with 1.6 Mq or the radius i? max 
of the maximum-mass Tolman-Oppcnheimer-Volkoff so- 
lution to characterize an EoS results in diagrams simi- 
lar to the upper left panel of Fig. [3J However, no clear 
trend can be found for the ejecta mass as a function of 
the maximum mass of nonrotating NSs. We therefore 
conclude that the NS compactness is the crucial EoS pa- 
rameter determining the ejecta mass. Indications of such 
a behavior were already observed in simulations for four 
simpli fied EoSs with an appro ximate temperature treat- 
ment (jHot okczaka et al.ll2013D . 

The dynamics of the merger explain why small NS radii 
lead to higher ejecta masses. For smaller -R1.35 the inspi- 
ral phase lasts longer and the stars reach higher velocities 
before they collide. The relation between the impact ve- 
locity and the NS radius is clearly seen in the lower left 
panel of Fig. [31 which displays the sum of the maxima 
of the coordinate velocities of the mass centers of the 
two binary components. The maximum of the coordi- 
nate velocity is reached shortly after the first contact, 
before the cores of the NSs arc decelerated by the col- 
lision. The clash of more compact NSs is more violent 
and more material is squeezed out from the collision in- 
terface, for which reason the negative correlation of the 
velocities with -R1.35 is reflected by a similar negative 
correlation of M C j CCta and -R1.35. Moreover, for smaller 
-R1.35 the central remnant consisting of the double cores 
rotates faster and the bounce and rebounce are stronger, 
i.e. the surface of the remnant moves faster and pushes 
away matter more efficiently. 

3.3. Influence of the approximate treatment of thermal 

effects 

A number of simulations of our survey (black and blue 
symbols) as w ell as calculations by other g roups (e.g. 
iRoberts et all 120111: iHotokezaka et all 120131) rely on an 
approximate description of thermal effects in the EoS, 
which requires the specification of an effective thermal 
ideal-gas index (see Sect. [5]). As can be seen in the up- 
per left panel of Fig. [31 the choice of the value for this 
ideal-gas index has a considerable impact on the ejecta 
mass; the simulations with T t h = 1-5 (blue symbols) yield 
generally more unbound matter. The reason is the re- 
duced thermal pressure support, which means that the 
two dense cores can approach each other more closely 



during the collision, which results in a more violent im- 
pact and shearing motion and thus in more material be- 
ing squeezed out from the collision interface and in a 
more powerful oscillation of the central remnant. This 
can be clearly seen by following the centers of mass of 
the two cores or the evolution of the central lapse func- 
tion. 

The optimal choice of T t h is a priori unclear and may be 
different for different EoSs. To address this issue we per- 
formed additional simulations (not shown in Fig. [3]) for 
temperature-dependent EoSs (SFHO, DD2, TM1, NL3) 
after reducing them to the zero-temperature sector (with 
the constraint of neutrino- less beta-equilibrium) and sup- 
plementing them with the approximate description of 
thermal effects using r t h = 1.5, 1.8 and 2. The compari- 
son with the fully consistent simulations reveals that gen- 
erally a choice of T t h = 1.5 yields the best quantitative 
agreement with only a slight underestimation of about 
10 per cent (20 per cent for the SFHO EoS), whereas the 
ejecta masses with T t h = 1.8 or r t h = 2 arc significantly 
too low compared to the fully consistent models (for the 
tested EoSs between 20 to 40 per cent for T t h = 2 and 
between 15 to 35 per cent for r t h = 1-8). 

The fact that a relatively low r t h reproduces the ejecta 
properties best contrasts the finding that a higher T t h (in 
the range between 1.5 and 2) has turned out to be more 
suitable for describing gravitational-wa ve features and 
the po st-merger collapse behavior (see IBauswein et al.l 
I2010aj) . i.e. the bulk mass motion of the colliding stars. 
The reason for this discrepancy is the density dependence 
of r t h, which drops from about 2 at supranuclear densi- 
ties to a bout 4/3 for den s ities be low re 10 11 g/cm 3 (see 
Fig. 2 in IBauswein et ail (|2010aD ). While the dynamics 
of the bulk mass of the merging objects, which is re- 
sponsible for the gravitational-wave production, is fairly 
well captured with a choice of r t h ~ 2, unbound fluid 
elements originating from the inner crust, where most of 
the ejecta stem from, encounter different density regimes. 
Consequently, the ejecta behavior cannot be well mod- 
elled with the r t h that is appropriate for high-density 
matter. We found the best compromise to be T t h ~ 1.5, 
but we stress that the results of simulations using an ap- 
proximate treatment of thermal effects should be taken 
with caution and do not need to be quantitatively reli- 
able in all aspects. 

3.4. Asymmetric binaries 

Even though most NS binaries are expected to be 
nearly symmetric systems with a total mass of about 
2.7 Mq, we investigate the mass ejection of an asymmet- 
ric setup to check whether also in this case an EoS de- 
pendence exists. The upper right panel of Fig. [3] displays 
the ejecta masses for simulations of asymmetric binaries 
with a 1.2 Mq NS and a 1.5 Mq NS. Again, the radius 
R1.35 of a nonrotating NS is used to characterize different 
EoSs. Here we restrict ourselves to EoSs which provide 
the full temperature dependence. In comparison to the 
symmetric binary mergers the amount of unbound mate- 
rial is significantly larger. The ejecta masses are about a 
factor of two higher than for the symmetric binaries with 
the same total binary mass. Also for asymmetric bina- 
ries a decrease of M e j with bigger -R1.35 is visible, but the 
scatter between models with similar -R1.35 is larger. The 
lower right panel of Fig. [3] shows the sum of the max- 
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Fig. 3. — Amount of unbound material for 1.35-1.35 Mq mergers (top left) and 1.2-1.5 Mq mergers (top right) for different EoSs 
characterized by the corresponding radius Hi, 35 of a nonrotating NS. Red crosses denote EoSs which include thermal effects consistently, 
while black (blue) symbols indicate zero-temperature EoSs that arc supplemented by a thermal ideal-gas comp onent with l\h = 2 (r t h = 1.5) 
(see main text). Small symbols represent EoSs which are incompatible with current NS mass measurements fDcmorcst ct al. 201Cj). Circles 
display EoSs which lead to the prompt collapse to a black hole. The lower panels display the sum of the maxima of the coordinate velocities 
of the mass centers of the two binary components as a function of -R1.35 for symmetric (bottom left) and asymmetric (bottom right) 
binaries. 



ima of the coordinate velocities of the mass centers of 
the two asymmetric binary components. As in the sym- 
metric case the two stars collide with a higher impact 
velocity if the initial radii of the NSs are smaller. 

Due to the asymmetry the dynamics of the merger pro- 
ceeds differently from the symmetric case (see Fig. HJ). 
Prior to the merging the less massive binary component 
is deformed to a drop-like structure with the cusp point- 
ing to the 1.5 Mq NS (top panels). After the stars begin 
to touch each other, the lighter companion is stretched 
and a massive tidal tail forms (middle left panel). The 
deformed 1.2 Mq component is wound around the more 
massive companion (middle panels). Also in the case of 
asymmetric mergers the majority of the ejecta originates 
from the contact interface of the collision, i.e. from the 
cusp of the "tear drop" and from the equatorial surface 
of the more massive companion, where the impact ab- 
lates matter (see top panels). Some matter at the tip 
of the cusp directly fulfills the ejecta criterion (top right 
panel), while the majority obtains an additional push 
by the interaction with the asymmetric, mass-shedding 



central remnant and the developing spiral arms (middle 
right and bottom panels). A smaller amount of ejecta of 
roughly 25 per cent originates from the outer end of the 
primary tidal tail (particles in the lower part of the top 
right panel). A part of this matter becomes unbound by 
tidal forces (at the tip of the tidal tail in the middle left 
panel) and the other fraction by an interaction with the 
central remnant (middle left panel). 

Figure [5] displays the distribution of the ejecta in a 
plane perpendicular to the binary orbit for the symmetric 
merger (left panel) compared to the asymmetric merger 
(right panel) for the last timesteps shown in Fig. [2] and 
Fig-IH respectively. A considerable fraction of the ejected 
matter is expelled with large direction angles relative to 
the orbital plane. For a timestep about 5 ms later the 
ejecta geometry is visualized (azimuthally averaged) in 
Fig. [6] excluding the bound matter. For both mergers 
the outflows exhibit a (torus or donut-like) anisotropy 
with an axis ratio of about 2:3. The velocity fields also 
show a slight dependence on the direction. 
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Fig. 4. — Same as Fig. [2] but for the asymmetric 1.2-1.5 Mq binary. Here the propagation speed of every 20th particle is indicated by 
an arrow and the side lengths of the panels differ from those of Fig. [2] In the upper panels the lower-mass star is identified by the black 
particles. 
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Fig. 5. — Distribution of the conserved rest-mass density (color-coded and logarithmically plotted in g/cm 3 ) in a plane perpendicular to 
the orbital plane for the 1.35-1.35 Mq merger (left) and the 1.2-1.5 Mq merger (right) with the DD2 EoS. The three-dimensional positions 
of unbound particles are projected into the cross-sectional plane, green and white referring to the origin from the one and the other NS. 
The velocity of every 20th particle is indicated by an arrwo. The time of the snapshots is given below the color bar of each panel. 




x[km] 20.0422 ms x [km] 19.9579 ms 

Fig. 6. — Ejecta geometry visualized by the rest-mass density (color-coded and logarithmically plotted in g/cm 3 ) excluding matter of 
the bound central remnant for the 1.35-1.35 Mq merger (left) and the 1.2-1.5 Mq merger (right) with the DD2 EoS. Density contours are 
obtained by azimuthal averaging. Arrows represent the coordinate velocity field where an arrow length of 200 km corresponds to the speed 
of light. The time of the snapshots is given below the color bar of each panel. 
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3.5. Binary parameter dependence 

The exploration of the full space of possible binary 
parameters is interesting for the determination of the 
highest and lowest possible ejecta mass for a given EoS 
and to understand the influence of the binary setup 
on the ejecta production. Such an investigation has 
been conducted o nly for one EqS (Shen) by Newto- 
nian calculations (iKorobkin et al.1 [2012; Rosswo g et al.1 
I20H iRosswod I2012D and within a relativistic frame- 
work (|Oechslin et al.l I2007D . which revealed quantita- 
tive differences between both approaches. Other sur- 
veys using different EoSs have been restricted to a lim- 
ited variation of the binary masses (iRoberts et al.lf20Tll 
iGorielv et all l20TTt iHotokezaka et al.l 120131 ) ~ Here we 
present the dependence of the ejecta mass on the mass 
ratio q = M1/M2 and the total binary mass M to t = 
M1 + M2 for a subset of EoSs employed in our study. The 
NL3, DD2 and SFHO EoSs are chosen because they are 
representative for the full set of possible EoSs: While the 
NL3 EoS is relatively stiff, resulting in -R1.35 = 14.75 km, 
the soft SFHO EoS produces rather compact NSs with 
-R1.35 = 11-74 km, and the DD2 represents an interme- 
diate case with R1.35 = 13.21 km. Figure [7] displays the 
amount of unbound matter as a function of the mass 
ratio q and total binary mass M to t for these three EoSs. 
The simulated binary configurations that form a differen- 
tially rotating NS are indicated in the figures by crosses, 
whereas systems leading to prompt black-hole formation 
(within about one millisecond after the first contact) are 
marked with circles (the ejecta properties are extracted 
10 ms after the merging). 

All three EoSs show qualitatively the same behavior. 
A clear trend of increasing ejecta masses with larger bi- 
nary asymmetry is visible. For symmetric binaries that 
do not undergo gravitational collapse there is also a slight 
tendency of higher total binary masses leading to more 
unbound material. The increase with M to t is more pro- 
nounced for asymmetric systems. The occurrence of a 
prompt collapse results in a significant drop of the ejecta 
mass. This is an important qualitative difference to New- 
tonian calculations, which cannot determine and follow 
the relativistic gravitational collapse. The threshold for 
the prompt collapse depends sensitively on the EoS, and 
soft EoSs lead to a collapse for relatively small M to t- 

For all EoSs the maximum ejecta masses (in all cases 
slightly below 0.02 M Q ) arc four to ten times higher 
than the amount of unbound matter of the symmetric 
1.35-1.35 Mq binaries. Here, the absolute differences be- 
tween the maximum and the minimum ejected mass for 
non-collapsing cases are larger for stiff EoSs like the NL3 
and less pronounced for soft EoSs like the SFHO. Soft 
EoSs yield steeper gradients of the ejecta mass in the 
binary parameter space, i.e. a certain variation in the 
binary parameters leads to a larger change of the ejecta 
masses than it is the case for a stiff EoS. To a good ap- 
proximation and ignoring cases with a prompt black hole 
formation, the setup with two stars of about 1.35 M is 
the system that produces the smallest amount of ejecta 
for the majority of investigated EoSs. The SFHO EoS 
is one of the exceptions, for which, e.g., the 1.2-1.2 M Q 
binary yields a factor two to three less ejecta than the 
1.35-1.35 Mq setup. 

Based on Newtonian calculations a fit formula for the 




Fig. 7. — Ejecta mass in A/q as function of the mass ratio q 
and the total binary mass Mtot for the soft SFHO EoS (upper 
panel), the intermediate DD2 EoS (center panel) and the stiff NL3 
EoS (bottom panel). In all panels the simulated binary setups arc 
marked by symbols. Crosses indicate the formation of a differ- 
entially rotating NS remnant, while circles identify configurations 
which lead to a direct gravitational collapse. 
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ejecta mass as fraction of M to t was pro posed as a func- 
tion o f r] = 1 - AMiM 2 /(M i + M 2 ) 2 in iKorobkin et al.l 
(|2012t) and iRosswod (|2012t) . Reviewing our data (even 
without the prompt collapse cases) we find a more com- 
plicated behavior and we can neither confirm the validity 
of the suggested fit formula nor find a generalization of 
it. This is not unexpected in view of the quantitative 
and qualitative differences between Newtonian and rela- 
tivists simulations discussed above. 

3.6. Folding with binary populations 

The dependence of the ejecta mass on the binary pa- 
rameters is essential to determine the total amount of 
ejecta produced by the binary population within a cer- 
tain time and thus to estimate the average amount of 
ejecta per merger event. The properties of the NS binary 
population are provided by theoretical binary evolution 
models, which still contain considerable uncertainties in 
many complexities of single star evolution and binary in- 
teract ion. Using the standard model of iDominik et all 
(|2012f) the folding of our results with the binary popu- 
lation yields an average ejecta mass per merger event of 
about 3.6 x 1CT 3 M Q for the NL3 EoS, 3.2 x 1CT 3 M 
for the DD2 EoS, and 4.3 x 1CT 3 M Q for the SFHO EoS. 
Therefore, the ejecta masses of the 1.35-1.35 Mq binary 
mergers give numbers for the three cases which approxi- 
mate the average amount of ejecta per merger event quite 
well (within 70 per cent for NL3, 3 per cent for DD2, 11 
per cent for SFHO). This finding is simply a consequence 
of the fact that the binary distribution is strongly peaked 
around nearly symmetric systems with M to t ~ 2.5 Mq 
so that the average ejecta mass is not sensitive to the 
larger ejecta production of asymmetric systems in the 
suppressed wings of the binary distribution. 

4. NUCLEOSYNTHESIS 
4.1. R-process abundances 

The potential of NS mergers to produce heavy r- 
process elements in their ejecta has been manifested 
by seve ral studies based on h y drodynamical simu- 
lations (iFreiburghaus et al.1 Il999t iMetzeer et all 120101: 
Roberts et all 1201 It iGorielv et all 120111: IKorobkin et all 
20121) . These investigations have consider ed only a few 
high-d ensity EoSs (two EoSs were used in IGorielv et all 
(|2011l )). Since the NS EoS affects sensitively the dynam- 
ics of NS mergers and thus the properties of the ejecta 
(amount, expansion velocity, electron fraction, tempera- 
ture) , we explore here the influence of the NS EoS on the 
r-process nucleosynthesis in a systematic way. 

For a selected, representative set of EoSs we extract the 
thermodynamical histories of fluid elements which get 
gravitational!/ unbound. For these trajec tories nuclear 
netwo rk calculations were performed as in IGorielv et all 
(|2011[) . where details on the reaction network, the tem- 
perature postprocessing and the density extrapolation 
beyond the end of the hydrodynamical simulations can 
be found. The reaction network includes all 5000 species 
from protons up to Z=110 lying between the valley of (3- 
stability and the neutron-drip line. All fusion reactions 
on light elements, as well as radiative neutron captures, 
photodisintcgrations, /3-decays and fission processes are 
included. The corresponding rates are based on experi- 
mental data whenever available or on theoretical prcdic- 
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Fig. 8. — Nuclear abundance pattern for the 1.35-1.35 Mq merg- 
ers with the NL3 (blue), DD2 (red) and SFHO (green) EoSs com- 
pared to the solar r-process abundance distribution (black). 
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Fig. 9. — Nuclear abundance pattern for the 1.2-1.5 Mq mergers 
with the NL3 (blue), DD2 (red) and SFHO (green) EoSs compared 
to the solar r-process abundance distribution (black). 

tions otherwise, as prescribed in th e BRUSLIB nuclear 
astrophysics library (|Xu et al. 2013) 

Figure[5]shows the final nuclear abundance patterns for 
the 1.35-1.35 M Q mergers described by the NL3 (blue), 
DD2 (red) and SFHO (green) EoSs. For every model 
about 200 trajectories were processed, which roughly cor- 
respond to about one tenth of the total ejecta. Compar- 
ing the final abundance distributions of the DD2 EoS for 
about 200 and the full set of 1000 fluid-clement histories 
reveals a very good quantitative agreement, which proves 
that a properly chosen sample of about 200 trajectories 
is sufficient to be representative for the total amount of 
unbound matter. 

The scaled abundance patterns displayed in Fig. [5] 
match closely the solar r-process composition above mass 
number A ps 140. In particular the third r-process peak 
around A = 195 is robustly reproduced by all models. 
Above mass number A ss 100 the results for the different 
NS EoSs hardly differ. For all three displayed models 
the peak around A 140 is produced by fission recy- 
cling, which occurs when the nuclear flow reaches fis- 
sioning nuclei around 280 No at the end of the neutron 
irradiation during the /3-decay cascade. The exact shape 
and location of this peak are therefore strongly affected 
by the theoretical modeling of the fission processes (in- 
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Fig. 10. — Nuclear abundance pattern for the 1.35-1.35 Mq merg- 
ers with the BSk20 (red) and BSk21 (blue) EoSs compared to the 
solar r-process abundance distribution (black). 

eluding in particular the fission fragment distribution of 
the fissioni ng nuclei) which are still subject to large un- 
certainties (jGorielv et al.l I2009T ) . Hence, the deviations 
from the solar abundance pattern between A 130 and 
A 170 are not unexpected, while the third r-process 
peak around A = 195 is a consequence of the closed neu- 
tron shell at N = 126, which is robustly predicted by 
theoretical models. Very similar results were obtained 
for NS merg er models per f ormed with the LS220 and 
Shen EoSs in lGorielv etall <|2011f ). 

In Fig.[9]the normalized abundance patterns are shown 
for asymmetric 1.2-1.5 Mq mergers employing the same 
representative EoSs as in Fig. [HI Again a very good agree- 
ment between the solar r-process abundances and the cal- 
culated element distributions above A w 130 is found for 
all three high-density EoSs. This confirms earlier find- 
ings that the binary mass ratio h as a negligible effec t on 
the abundance yield distribution (|Gorielv et al.ll20il"l) . It 
also confirms that the ejected abundance distribution is 
rather insensitive to the adopted EoS. 

Besides the three temperature-dependent EoSs consid- 
ered above we conducted network calculations also for 
merger models computed with zero-temperature EoSs 
supplemented by an approximate treatment of thermal 
effects with a I\h = 2 ideal-gas component (BSk20, 
BSk21). In this case the tem perature is estimated follow- 
ing a procedure described in lEtienne et al.1 (|2008| ). which 
converts the specific internal energy to temperature val- 
ues. Doing so it is assumed that the energy of the thermal 
ideal-gas component is composed of the thermal energy 
of an ideal nucleon gas and a contribution from ultrarel- 
ativistic particles (photon, possibly electrons, positrons 
and neutrinos). 

The network calculations for the BSk20 and BSk21 
EoSs yield an abundance pattern above A rts 130 very 
similar to the other fully temperature-dependent EoSs 
(see Fig. [TO]) . Differences between the fully consistent 
models and the simulations with approximate tempera- 
ture treatment are found below mass number A « 50, 
where the calculations with the BSk EoSs yield a lower 
amount of elements with 5 < A < 50 but a higher mass 
fraction of hydrogen, deuterium and helium. The reason 
is the higher temperatures found with the BSk EoSs at 
the beginning of the network calculations which lead to 



a reduced recombination of nuclcons and a-particles and 
consequently a smaller production of heavier nuclei. In 
this respect the conditions in the outflows of these mod- 
els resemble the situation in the neutrino-driven winds 
of core-collapse supcrnovae but for significantly higher 
neutron excesses. 

Overall, it is reassuring that the BSk models yield 
a similar abundance pattern of r-process elements, al- 
though these calculations rely on an approximate incor- 
poration of thermal effects and a rough estimate of tem- 
peratures. It is an important finding of our work that 
r-process elements are robustly produced for a represen- 
tative, diverse sample of high-density EoSs and that the 
outcome is insensitive to the exact initial temperature 
conditions and the binary setup. 

4.2. Merger rates 

The above-mentioned variations in the production of 
light elements are also reflected in the fraction of the 
ejecta which end up as r-process elements. In Tab. Q] 
a clear difference is observed for temperature-dependent 
EoSs (NL3, DD2, SFHO) and the BSk EoSs with ap- 
proximate temperature treatment. While in the former 
cases about 96 to 99 per cent of the ejecta are converted 
to r-process elements, only 93 to 95 per cent of the ejecta 
are processed to r-process material in the latter models. 
This variation is a consequence of the temperature his- 
tory affecting the production of light nuclei, as discussed 
above. 

Despite these differences, it is justified to assume that 
almost the total amount of ejecta is converted to heavy 
r-process elements, also considering the uncertainties in 
the determination of the exact ejected masses from simu- 
lations. Furthermore, in Sect.[3]it was argued that the to- 
tal amount of ejecta, and thus the r-process material pro- 
duced per event by the population of NS binaries, can be 
well represented by the yield of the 1.35-1.35 Mq merger. 
This allows an important consistency check by compar- 
ing the theoretically expected production to the observed 
amount of r-process matter with A > 130 in the Galaxy , 
which is estimated to be about 4 x 10 3 M Q ()Qianll2000l) . 
In order to produce this amount of heavy r-process ele- 
ments within the Galactic history for about 10 10 yr, one 
requires a merger rate of 4 x 10 /yr if every coalescence 
ejects on average 10 -3 Mq, which in our EoS survey 
corresponds to the lower bound on the ejecta mass of 
1.35-1.35 Mq mergers (see upper left panel of Fig. |3]). 
Similarly, assuming that NS mergers are the dominant 
source of heavy r-process elements, the upper bound of 
10~ 2 Mq for the ejected mass from 1.35-1.35 Mq binaries 
(i.e for EoSs with small -R1.35) would be compatible with 
a merger rate of4x 10 _5 /yr. These rate estimates lie in 
the ballpark of th eoretical predictions ranging from 10 -6 
to 10 -3 per year (jAbadie et al.ll2010l) . This implies that 
all EoSs of our survey are compatible with NS mergers 
being the dominant or a major source of r-process ele- 
ments. 

This conclusion on the merger rate may be tested 
against future observations, in particular by multiple 
gravitational-wave detections or frequent observations 
of electromagnetic counterparts. Our work emphasizes 
that, in addition to a more accurate merger rate, infor- 
mation on the high-density EoS is needed to shed light 
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on NS mergers as a major source of r-process elements. 
More specifically, a merger rate of about 4 x 10 _5 /yr may 
imply cither that nearly all heavy r-process elements arc 
made by NS mergers in the case of a soft high-density 
EoS with small NS radius -R1.35, or that only a tenth of 
the observed r-process material originates from NS bina- 
ries if a stiff NS EoS with large -R1.35 is confirmed. 

Inversely, considering the robustness of the r-process 
nucleosynthesis in NS mergers, one can infer that, assum- 
ing a minimal production of ~ 10 -3 M of r-nuclei per 
event a constant merger rate during the life of the Galaxy, 
this rate cannot be higher than roughly 4 x 10 _4 /yr; oth- 
erwise the galactic r-process material would be more than 
presently observed in the Galaxy. This b ound is com- 
parab le to the "optimistic" limit given in lAbadie et al.l 
(|2010T l . Thus, the r-process element content in the 
Galaxy establishes further, independent evidence for an 
upper limit on the merger rate below ~ 4 x 10 -4 /yr. 
A restrictio n to soft EoSs , e.g. from other physical con- 
straints like lHebeler et all (|2010f ). would lead to a smaller 
upper limit for the event rate of NS mergers. 

4.3. Actinide production ratios and stellar chronometry 

Some of the heaviest long-lived radioactive nuclei 
produced by the r-process can be used as nucleo- 
cosmochronometers. In particular the abundance ra- 
tios of thorium to europium and thorium to uranium 
have been proposed for estimating the age of the old- 
est stars in our Galaxy. More specifically, a sim- 
ple comparison of the observed abundance ratio with 
the production ratio can provide an age estimate of 
the c o ntaminated object ( Butcherl 119871: lFrancois et all 
1991 ICowan et all 1199a IGorielv fc Clerbauxl 1199a 



Gorielv fe Arnouldll200lf[Cavrel et al.ll2001l: iFrebel et all 
20071 : iSneden et all 120081) . In addition, if we consider 
low-mctallicity stars polluted by a small number of nu- 
cleosynthetic events that took place just before the for- 
mation of the stars, the age of the star can be estimated 
without calling for a complex model of the chemical evo- 
lution of the Galaxy. The major difficulty of the method- 
ology is therefore related to the theoretical estimate of 
the r-production ratio and the corresponding uncertain- 
ties of astrophysics and nuclear physics origin that may 
affect this prediction. 

In this respect, the 232 Th to 238 U chronometry has 
been shown to be relatively robust, in particular in com- 
parison with the Th/Eu chronometry, i.c to be less af- 
fected by the still large astrophysics and nuclear physics 
uncertainties aff ecting our understanding of the r-process 
nucleosynthesis (fGoriely fe Clerbauxlll99a ICavrel et all 
[20011: IFrebel et alll2007l) . In Tab.fTlwe provide the pro- 
duction ratios of the 232 Th to 238 U isotopes based on 
our NS merger simulations. Assuming that r-process 
enhanced metal-poor stars were enriched by one or a 
few NS merger events, we can derive ages within this 
scenario. From the observed ratio of log (U/Th) obs = 
-0.9 4 ± 0.09 for the metal-po or star CS31082- 
001 ilCavrel et al.l l200l (see also IGorielv fc Arnouldl 
((2001 ) for updated values), we compute the age as At = 
21.8 [log (U/Th) - log (U/Th) ob J = 15.7 Gyr with the 
production ratio log (U/Th) = -0.22±0.01 (see Tab.H]). 
The observational uncertainty of 0.09 dex in this case 
dominates the error on the age estimate since it amounts 
to about 2.0 Gyr, while the theoretical uncertainties as- 



sociated with the different EoSs (Tab.[T]) give an 0.2 Gyr 
error only. Significant additional uncertainties stem from 
the nuclear physics aspects of the r-process nucleosynthe- 
sis; those will be studied within the NS merger model in 
a forthcoming paper. 

The derived age of this h alo star lies within the ball- 
park of other age estimates ([Cavrel et al.ll200"lT ). For the 
metal- poor star HE 1523 -0901 (log (U/Th) Q&s = -0.86 ± 
0.13) (|Frebel et al.|[2007h our 232 Th to 238 U production 
ratio implies an age of about 14.0 ± 2.8 Gyr, which is 
also w ithin the range of other calculations (|Frebel et all 
120071) . 

Tab. Q] also lists the production ratios of 232 Th to 
235 U as well as Th to Eu. For age estimates the 
Th/Eu chronometer has been widely used, although 
it remains highly sensitive to all types of uncertain- 
ties. In this case, the stellar age is derived from 
At = 46.7 Gyr [log (Th/Eu) - log (Th/Eu) ob J, so that 
a 25% error on the production or observed Th/Eu ra- 
tio gives rise to an uncertainty of about 5 Gyr on 
the stellar age. Special care of the associated uncer- 
tainties should the refore be taken when ap plying this 
chronometer pair (IGorielv fc Arnouldl l200l . Consid- 
ering our production ratio of about log (Th/Eu) obs = 
-0.20 ±0.06 (see Tab. [J), we find an age of At = 17.7 ± 
2.8 Gyr for HE 1523-0901 (log (Th/Eu) obs = -0.58, 
with an additional 4.8 Gy r uncertainty based on obser- 
vation (IFrebel et al.l 120071) ). We stress again here that 
as long as the r-process site remains unidentified, cor- 
responding uncertainties of the production ratios have 
to be taken into account as well as uncertainties arising 
from the nuclear physics input in network calculations. 
We also point out the possibility that there is a mea- 
surable event-to-event variation in the production ratios, 
for instance in the case of NS mergers caused by the 
unknown binary configuration (see Tab. [T]), but also for 
other sites a progenitor dependence cannot be excluded. 

The production ratios summarized in Tab.[T]are of par- 
ticular importance since NS mergers may well be a major 
source of r-process elements and current supernova mod- 
els cannot provide suitable conditi ons for the formation 
of the heaviest r-process elements dHoffman et all 120081: 
Janka et all 120081: iRoberts et al.l 120101 : iHudepohl et al l 



20101: iFischer et alli2oTot iWanaio et aLll201lD . The rela 



tively reliable age estimates from the 232 Th to 238 U ratio 
are compatible with the age of the universe, and thus NS 
mergers cannot be excluded as the source of the contam- 
ination of the considered metal-poor stars. 

4.4. Nuclear heating 

Another outcome of our nucleosynthesis calculations 
is the determination of the heating due to radioac- 
tive decays in the ejecta. This is in particular impor- 
tant because the radioactive heating provides the en- 
ergy source for an optical coun t erpar t associated with 
NS mergers (see ILi fc Paczvnskil (fl998h . iKulWl ([2005) 
and IMetzger et all (|2010f ) and Sect. ®. Our calcula- 
tions allow us to check the robustness and general be- 
havior of the heating rate. In all models the heating 
rate due to radioactive decays (beta-deca ys, fission and 
alpha- decays) looks similar (see Fig. 3 in IGorielv et all 
(|201lD ). For instance at the time t pea k, when the lu- 
minosity of the optical transient reaches its maximum 
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TABLE 1 

Nucleosynthesis calculations 



Model Mx-M 2 


M ej 


^pcak 


■A^r— process / -^ej 


232^/238^ 


232rp n y235jJ 


Th/Eu 




/ 




(io- 3 M ) 


(d) 










(MeV/A) 




NL3 1.35-1.35 


2.09 


0.171 


0.989 


1.644 


1.072 


0.695 


3.34 


1.8 x 10" 6 


DD2 1.35-1.35 


3.07 


0.189 


0.980 


1.671 


1.080 


0.627 


3.13 


1.6 x 10~ 6 


SFHO 1.35-1.35 


4.83 


0.228 


0.991 


1.642 


1.039 


0.579 


3.17 


1.4 x 10" 6 


NL3 1.2-1.5 


7.95 


0.338 


0.964 


1.670 


1.112 


0.714 


3.36 


1.3 x 10~ 6 


DD2 1.2-1.5 


8.79 


0.354 


0.986 


1.697 


1.109 


0.658 


3.11 


1.2 x 10~ 6 


SFHO 1.2-1.5 


13.39 


0.418 


0.974 


1.685 


1.085 


0.543 


3.12 


1.1 x 10~ 6 


BSk21 1.35-1.35 


3.36 


0.162 


0.948 


1.660 


1.023 


0.704 


2.97 


1.6 x 10~ 6 


BSk20 1.35-1.35 


4.68 


0.195 


0.931 


1.689 


1.021 


0.698 


2.96 


1.5 x 10" 6 



Note. — Selected models for which nucleosynthesis calculations were performed. M e j is the amount of unbound matter, whereas t pea k 
is the peak time of an optical transient associated with a NS merger (see Sect. |5j. The fourth column gives the fraction of the ejecta 
which is processed into r-process elements. The production ratios of certain elements and isotopes are provided in the fifth to seventh 
columns. £?h oa t denotes the total amount of energy released by radioactive decays (without neutrino energy). The factor / approximates 
the radioactive heat generation around the time of the optical peak luminosity relative to the rest-mass energy of the ejecta (see text). 



(typically several hours; see Sect. [5]), the heating rate 
varies from 3 x 10 10 erg/g/s to 1 x 10 11 erg/g/s for the 
cases where detailed nucleosynthesis calculations were 
made. This implies that the heating efficiency / = 

Q(ipeak)ipeak/M3jC 2 that enters the es timates of opti- 
cal em issio n properties (see Sect . [5] and iLi fe Paczvnskil 
dHH) and IMetzger et all (|2010D V varies between 1.1 x 
10 -6 and 1.8 x 10~ 6 (see Tab. [J). One observes a moder- 
ate dependence of / on the EoS and the mass ratio. This 
is caused by the longer duration i pea k of the emission 
peak for larger ejecta masses in models with soft EoSs or 
asymmetric binaries. Overall / « 1.5±0.3 x 10 -6 seems 
a fair approximat i on, w hich is half of the value suggested 
in IMetzger et al.1 (|2010l ) for a longer peak time of about 
one day. As detailed in Sect. [5J we expect from our rel- 
ativistic merger simulations shorter peak times in the 
range of 2 to 7 hours for symmetric binaries. The total 
amount of energy released by radioactive decays (with- 
out neutrinos) is in the range of about 3.2 ± 0.2 MeV per 
nucleon (see Tab.[T|). 

5. OPTICAL COUNTERPARTS AND RADIO REMNANTS 

The radioactive decay of the synthesized r-process ele- 
ments generates heat, which is d eposited in the ejecta 
and powers an optical display (ILi fe Paczvriski| 



Kulkarnill200i| IMetzger et al.ll2010t iRoberts et al.l 
Gorielv et all 120111) . This electromagnetic counterpart 



1998 



2011 



of a NS merger is potentially observable with existing 
and upcoming optical surveys such as the Palomar Tran- 
sient Factory, the Synoptic All Sky InfraRed Survey, 
the Panoramic Survey Telescope and Rapid Response 
Syst em, and the Large Synopti c Sur vey Telescope (see 
e.g. iStrubbe fe Quataertl (|2009f) and IMetzger fe Bergerl 
(|2012f ) for a compilation o f certain characteri stics of these 
facilities). We also refer to lKann et al.l (j 2 1 If ) for a report 
on attempts to detect signatures of a radioactivcly pow- 
ered transient in the light curve following a gamma-ray 
burst. The detection of such optical signals would pro- 
vide valuable information on ejecta properties and the 
sky position of an event. As detailed in the next section, 
the peak luminosity, peak time, peak width, and effective 
temperature depend on the amount of ejecta and the ex- 
pansion velocity. From an observation one might there- 
fore derive these ejecta characteristics, which arc other- 
wise only accessible by numerical simulations. Detecting 
radioactively powered emission and determining ejecta 



masses will consolidate the role of NS mergers for the 
enrichment of the Galaxy with heavy r-process elements. 
A precise localization of an unambiguously identified op- 
tical transient of a NS merger event will help improving 
the sensitivity of gravitational-wave detections, and will 
provide information about the host galaxy or environ- 
ment. Moreover, if the distance scale of the events can 
be constrained, better observational limits on the merger 
rate will become available. 

5.1. Model 

The bolometric peak luminosity of an optical transient 
associated with a NS coalescence can be estimated by 



^pcak 



5xl0 41 erg/s 



_t_\ /jl 

10" 6 J VO.lc 
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with the average outflow velocity v, the ejecta mass 
M G j and the he ating efficien c y / already intro- 
duced in Sect PI (lArnettl 119821 : (Li fc Paczvnskl 119981 : 
IMetzger et al.M2010f ). 

The time of the peak luminosity and the effective tem- 
perature at the time of the maximum luminosity can 
also be expressed as functions of the outflow velocity, 
the ejecta mass, and /: 



ipcak ~ 0.5 d 



T peak « 1.4xl0 4 K 



v 



-1/2 



M. 



10- 2 Af, 



1/2 



(3) 



10" 



1/4 



0.1c 



-1/8 



10" 2 M n 



-1/8 



(4) 

(see ILi fc Paczvnskl (|1998l) and IMetzger et all (|2010D ). 
Within this model the width Ai poa k of the lumi- 
nosity peak is propor tional to i pca k (jArnettJ 119961 : 
iKasen fc Wooslevl I2009D . On the basis of the one-zone 



model in IGorielv et alT(|2^ll we find that the full width 
at half maximum can be very well approximated as 



2.5i pea k- 



(5) 



A£ pca k 

The above formulas can be understood by some general 
considerations. For larger M e j or smaller v the ejecta 
need a longer time to become transparent. If the out- 
flow gets optically thin at a later time, expansion cooling 
reduces the effective temperature at £ peak - Since, how- 
ever, the increase of the emission radius i? pca k ~ v x i peak 
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dominates the decrease of the temperature in the Stcfan- 
Boltzmann law, L pca k oc ^peak-^pcakJ the peak luminos- 
ity increases with larger M c j and higher v. The simple 
scaling laws of Eqs. [2] to 2] for the prope rties of optical 
transi e nts are confirm e d by c alcul ations inlMetzger et all 
((201?]) . IRoberts et all (l201lh and lGorielv et all (|2011D . 

5.2. EoS dependence 

The panels on the left side of Fig. [TT] display the peak 
luminosity, the peak time and the effective temperature 
of electromagnetic counterparts of 1.35-1.35 Mq merg- 
ers for different EoSs, which arc characterized by the 
radii of the corresponding 1.35 Mq NSs. A clear de- 
pendence of the optical display on the compactness of 
the NSs can be seen, with soft EoSs yielding brighter 
transients, which peak on longer timescales with a lower 
effective temperature. The relatively clear relations are 
mainly a consequence of the strong EoS impact on M e j 
(Fig. [3]), while the average outflow velocity varies only by 
a factor of 3 (see Fig. [12"]) . The average expansion veloc- 
ities show the tendency of being higher for EoSs which 
yield smaller -ffi.35 (see Fig. [12]). This is consistent with 
the reasoning used before that more compact NSs lead 
to more violent collisions. However, the relation between 
the average outflow velocity and -ffi.35 is not very tight. 
For symmetric binaries the outflow velocities (measured 
10 ms after the merging at which time the asymptotic 
values are fairly well determined) vary (with larger scat- 
ter) between 0.16 and 0.45 times the speed of light. 

Also for asymmetric 1.2-1.5 Mq binaries an EoS depen- 
dence of optical counterpart properties is observed (right 
side of Fig. [TTj) . As in the case of equal-mass mergers, 
EoSs with smaller -R1.35 lead to more mass ejection and 
therefore to more luminous events with longer i pea k and 
lower T e fj. Compared to symmetric binaries asymmetric 
setups generally produce brighter transients, which reach 
their peak luminosities on longer timescales and therefore 
lower effective temperatures, because their ejecta masses 
are higher whereas the expansion velocities are compa- 
rable with those of symmetric systems described by the 
same EoS (see Fig. [12]). For all considered quantities 
the asymmetric models exhibit a milder EoS dependence. 
Plotting the peak luminosity as a function of the binary 
mass ratio and the total binary mass for the NL3, DD2 
and SFHO EoSs reveals a qualitatively similar behavior 
as the ejecta masses shown in Fig. [7] 

The relations shown in Fig. [TT] suggest the possi- 
bility to constrain NS radii and thus the high-density 
EoSs from observations of optical transients associated 
with NS mergers. Optimally such a detection could 
be supplemented by a gravitational- wave measurement, 
which provides the involved binary masses, the distance, 
and the merger time. But even without an associated 
gravitational-wave signal the observable features of a 
transient may have the potential to yield constraints for 
the NS EoS. The combinations of L pe ak, ipeak, and T ff 
vary systematically with the NS properties. For instance 
a low peak luminosity, a small peak width, and a high 
effective temperature imply a large NS radius. 

5.3. Implications for observations 

Symmetric 1.35-1.35 Mq binaries are predicted to be 
the most common configurations and thus arc likely to 



be the ones first and most frequently observed. Unfor- 
tunately, the 1.35-1.35 Mq systems yield the smallest 
ejecta masses and thus the lowest luminosities and peak 
timescales. The peak widths are important to estimate 
the prospects of blind searches, whereas the peak time 
sets the scale for the response time after a gravitational- 
wave trigger. 

The nearly complete coverage of EoS possibilities by 
our survey allows us to determine the possible range 
of signal properties of optical transients associated with 
NS merger events. From our survey we find that the 
optical peak luminosity of a 1.35-1.35 Mq NS merger 
should be expected to be between about 3 x 10 41 erg/s 
and 14 x 10 41 erg/s corresponding to absolute bolometric 
magnitudes of M = -15.0 and M = -16.7. The peak 
times range from only 2 hours to 7 hours, and the dura- 
tion of the emission is expected to be between 4.8 hours 
and 18 hours depending on the high-density EoS. Note 
that the ballpark of our models yields fainter and shorter 
transients than typical estimates based on Newtonian 
models, which in general obtain hig her ejecta masses and 
lower average expans i on velocities (IRoberts et all 120111: 
Korobkin et all l2012t iRosswog et al.l 120121: iPiran et all 
2012tlRosswodl2012|) (see also the discussion in Sect. [3]). 
While for the peak luminosity these differences lead to 
partially compensating effects, the timescales are more 
strongly affected. For symmetric binaries even the max- 
imum peak time of about 7 hours found in our sample 
is well below most predictions based on Newtonian mod- 
els (iMetzger et alj|2010t IRoberts et al.ll2011t|Piran et al.l 
120121 : IRosswog et al.ll20lf ~ 

As shown in [6] during the early stages of the expan- 
sion the ejecta exhibit a fair asymmetry between polar 
and equatorial directions. The forumlas used in this sec- 
tion for estimating the properties of optical counterparts 
assume spherically symmetric outflows. It remains to be 
explored whether the donut-like shape visible in Fig. [6] 
persists at late times or whether the outflow becomes 
more symmetric when the peak of the optical display 
occurs. Multidimensional radiation transport calcula- 
tions coupled to long-term hydrodynamical simulations 
are required to determine to which extent the simplified 
emission model provides reliable estimates of the observ- 
able properties of the electromagnetic transients in de- 
pendence on the observer direction. 

5.4. Radio flares 

Another potentially observable phenomenon connected 
with NS mergers is radio emission that is produced by 
the interaction of the outflowing ma t erial with the am- 
bient medium dNakar fc Piranl [20 lit IPiran et all 120121: 
IRosswog et al.ll2012h . The ejecta properties which deter- 
mine the appearance of a radio remnant are the kinetic 
energy and the outflow velocity. The peak flux density 
was computed to be proportional to the total kinetic en- 
ergy Ek in of the outflow and to the (initial) outflo w ve- 
locity v to a power of about 2.5 (|Piran et al.ll2012l ). For 
symmetric mergers we find values for the kinetic energy 
between 6 x 10 erg and 10 51 erg (Fig. [T3]). whereas the 
average outflow velocities vary from 0.16 to 0.45 times 
the speed of light (Fig. IT2"]) . The kinetic energy scales 
well with the expansion velocity, i.e. the models with 
the highest outflow velocities yield also the highest ki- 
netic energies, and configurations with smaller v result 
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Fig. 11. — Estimated properties of the optical transients for symmetric (1.35-1.35 Mq) mergers (left panels) and asymmetric (1.2-1.5 Mq) 
mergers (right panels) for different EoSs characterized by the NS radius K1.35. The symbols have the same meanings as in Fig. [3] The 
top panels show the bolometric peak luminosity, the middle panels the corresponding peak timescale, and the bottom panels the effective 
temperature at the time of the peak luminosity. 
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Fig. 12. — Average ejecta expansion velocity for 1.35-1.35 Mq 
mergers (with symbols analog to Fig. [3j and for 1.2-1.5 Mq (red 
squares) for different EoSs characterized by the corresponding ra- 
dius -R1.35 of a nonrotating NS with a mass of 1.35 Mq. 
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Fig. 13. — Kinetic energy of the ejecta for 1.35-1.35 Mq mergers 
(with symbols analog to Fig. [3j and for 1.2-1.5 Mq (red squares) 
for different EoSs characterized by the corresponding radius R1.35 
of a nonrotating NS with a mass of 1.35 Mq. 

in lower kinetic energies. This implies that the peak flux 
of radio remnants is uncertain by a factor of 200 because 
of the variations in the kinetic energy and expansion ve- 
locity associated with the incomplete knowledge of the 
NS EoS. Similar values are found for asymmetric merg- 
ers. 

6. SUMMARY AND CONCLUSIONS 

We have performed relativistic hydrodynamical simu- 
lations of NS mergers to investigate the mass ejection, 
the nucleosynthesis outcome, and the properties of as- 
sociated optical transients. The main goal of this study 
was to explore the systematics of the EoS dependence of 
these aspects by employing a large set of candidate EoSs, 
while focussing mostly on binaries with two 1.35 Mq NSs 
and on asymmetric 1.2-1.5 Mq systems. 

The unbound ejecta mass is strongly affected by the 
adopted EoS. We find that the NS compactness is the 
crucial EoS parameter determining the ejecta properties. 



Using the radius -R1.35 of a nonrotating 1.35 Mq NS to 
characterize different EoSs we find that "softer" EoSs 
which yield more compact NSs tend to produce more 
ejecta. The ejecta masses are between 10 -3 Mq and 
1.5 x 10 -2 Mq depending on the EoS and binary mass 
ratio (Figs. |3] and [7]). Most of the unbound material 
originates from the contact interface between the collid- 
ing stars for symmetric as well as asymmetric binaries. 
In the latter case about 25 per cent of the mass ejec- 
tion are shed off the outer end of the spiral arm like tail 
into which the lower-mass component is stretched dur- 
ing its final approach to collision with the more massive 
companion. 

A qualitative and quantitative agreement of our SPH 
simulations (employing the conformal flatness approx- 
imation) with fully relativistic grid-based simulations 
is found, whereas considerable differences compared to 
Newtonian models concerning the origin and the amount 
of ejecta are observed. The pronounced spiral arms, for 
example, which form in Newtonian simulations during 
the merging of symmetric binaries and whose mass strip- 
ping dominates the ejecta, are absent in relativistic merg- 
ers of equal-mass NSs. Newtonian models therefore tend 
to produce considerably higher ejecta masses. 

When temperature effects are mimicked by adding a 
thermal ideal-gas component with a constant ideal-gas 
index T t h to EoS models which are provided as zero- 
temperature barotropes, the best match of the ejecta 
masses with fully consistent calculations is achieved for 
a relatively low value of l\h = 1-5. This is in con- 
flict with values of 1.8 or 2 which have been widely 
used and work well for gravitation al-wave determina- 
tions (|Bauswein et all [20121 [20l0al ). This can be un- 
derstood from the fact that the gravitational-wave sig- 
nal is produced by the bulk of the merger mass in the 
high-density regime, whereas the ejecta depend on the 
thermodynamics of lower-density matter expanding away 
from the colliding stars and being accelerated by pressure 
forces. 

The binary parameters have qualitatively the same in- 
fluence on the ejecta masses for all investigated EoSs. A 
larger binary mass asymmetry leads to a strong increase 
of the mass of unbound matter, whereas a higher total bi- 
nary mass results in larger ejecta masses for asymmetric 
mergers but only a weak increase of M e j for symmetric 
systems (Fig. [7]). The occurrence of a prompt collapse of 
the merger remnant is associated with a significant drop 
in the ejecta mass. For a given EoS (and no prompt 
collapse) the smallest amount of ejecta is to a good ap- 
proximation produced by 1.35-1.35 Mq binaries, while 
the ejecta mass can be up to a factor ten higher for other 
binary configurations. For soft EoSs the ejecta masses 
show steeper gradients in the binary parameter space 
(q, Mtot). 

Nuclear network calculations show that r-process ele- 
ments with mass numbers above 130 are robustly pro- 
duced in the ejecta of NS mergers for all investigated 
EoSs. The vast majority of ejected material is fission re- 
cycled producing a final mass-integrated abundance pat- 
tern that resembles closely the solar composition of r- 
process elements. The robustness with respect to varia- 
tions of the high-density EoS confirms that NS mergers 
are a very promising source of r-process elements. 

For some EoS models which do not provide the tern- 
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able features (luminosity, peak timescale, effective tem- 
perature) and NS radii we propose that optical obser- 
vations of transients associated with NS mergers could 
yield valuable constraints on the NS EoS. 

The very broad range of possibilities included in our 
EoS sample allows us to bracket the expected range of 
signal features of optical counterparts associated with 
NS mergers. Optical transients of 1.35-1.35 M Q merg- 
ers should (at least) reach an absolute bolometric peak 
magnitude between -15.0 and -16.7 (3 x 10 41 erg/s and 
14 x 10 41 crg/s). Depending on the high-density EoS 
the peak times vary from 2 to 7 hours, implying du- 
rations of about 4 to 18 hours, whereas effective tem- 
peratures between 1.3 x 10 4 K and 1.9 x 10 4 K can 
be expected. We emphasize that the peak luminosi- 
ties, peak times, and peak widths of the optical coun- 
terparts are found to be considerably lower in our anal- 
ysis compared to earlier investigations based on New- 
tonian models dMetzeer et alj|2010l iRoberts et all 1201 It 



peraturc dependence consistently, the hydrodynamical 
simulations are based on an approximate treatment of 
thermal effects and the nucleosynthesis calculations rely 
on an estimation of the temperature in the ejecta. Also 
in these cases r-process elements are produced with solar 
distribution, showing the insensitivity to the exact tem- 
perature value at the beginning of the network calcula- 
tions. Moreover, the abundance patterns are not affected 
by asymmetries in the binary setup. The amazing insen- 
sitivity of the outcome of the nucleosynthesis processes 
can be understood as a consequence of the large neutron 
excess in the ejected inner-crust material of the merging 
NSs, which allows for fission recycling in essentially all 
considered conditions. 

By folding with the binary population we identify the 
results of 1.35-1.35 Mq mergers as a good approxima- 
tion for the average ejecta mass per merger event. The 
main uncertainty in the average ejecta mass per merger 
is therefore associated with the incomplete knowledge 
of the NS EoS, which implies variations in the average 
ejecta mass of about a factor of 10. The observed abun- 
dance of r-process matter in the inventory of our Galaxy 
can be accounted for with a NS merger rate that is com- 
patible with current predictions based on population syn- 
thesis and pulsar observations (lAbadie et al.l [20Toh . Fi- 
nal conclusions, however, require a more precise determi- 
nation of the merger rate, e.g. by gravitational-wave de- 
tections or observations of electromagnetic counterparts. 
Our work implies that in addition to more accurate in- 
formation on the merger rate also better constraints on 
the high-density EoS are needed to decide whether NS 
mergers are a major (or the dominant) source of r-process 
elements. 

Moreover, our simulations in combination with esti- 
mates of the Galactic r-process material provide inde- 
pendent evidence that the Galactic merger rate cannot 
be higher than approximately 4 x 10 -4 events per year if 
NS mergers should not overproduce heavy r-nuclei com- 
pared to observations. 

The nucleosynthesis calculations of our survey also pro- 
vide important information on the production ratios of 
certain isotopes which are used for nuclcocosmochronom- 
etry. For instance the ratio of 232-thorium to 238- 
uranium is found to be about 1.65 with only small varia- 
tions depending on the high-density EoS and the binary 
configuration (Tab.[lJ. Using this result we derived ages 
of metal-poor stars which are consistent with other age 
estimates. This implies that NS mergers are not excluded 
as r-process element sources for metal-poor stars and the 
production ratios provided here for the first time in the 
NS merger context should be taken into account in stel- 
lar age estimates as long as mergers cannot be excluded 
as r-process sites in the early Galactic history. 

Just as the ejecta masses of binary NS mergers ex- 
hibit a strong sensitivity to the properties of the nuclear 
EoS and thus to the radius -R1.35 of the merging stars, 
we also predict the optical transients powered by the ra- 
dioactive energy release in the ejecta to depend on the 
compactness of the binary components. EoSs which lead 
to smaller NS radii produce more ejecta and therefore 
cause brighter optical counterparts, which peak on longer 
timescales with longer durations and with lower effective 
temperatures. On the basis of our extensive survey of 
EoSs, which suggests clear correlations between observ- 



iPiran et alJl2012tlRosswog et al.ll2012t) . The reduction is 
a consequence of the smaller ejecta masses especially for 
symmetric binaries. Because of the shorter peak time 
and duration of the optical transient suggested by rcl- 
ativistic merger results, we also find a smaller fraction 
/ of radioactive decay energy relative to the rest-mass 
energy of the ejecta. While Newtonian merger mod- 
els yi eld / ~ 3 x 10~ 6 a t the time of the luminosity 
peak (jMetzger et alJl2010f ) we obtain / - 1.5 x 10~ 6 with 
little sensitivity to the EoS and the binary parameters. 

For different EoSs considerable differences are also 
found in the properties of radio remnants. Based on our 
sample of models we estimate an uncertainty of up to a 
factor of 200 for the theoretical predictions of the bright- 
ness of these events. 

Future work should address a variety of issues. The hy- 
drodynamical models of NS mergers should include mag- 
netic fields and the effects of neutrino interactions should 
be explored. Our results also need to be confirmed by 
fully relativistic merger simulations. The properties of 
emitted electromagnetic radiation should be computed 
for detailed multi-dimensional outflow models including 
the corresponding nuclear network calculations to deter- 
mine the composition and heating. Radiative transfer 
calculations will have to be performed, employing appro- 
priate opacities of r-process elements, to study the obser- 
vational appearance of the potentially anisotropic ejecta 
dependent on the viewing direction. In this context it 
will also be important to determine the contribution of 
mass ejection from the secular evolution of the merger 
remnant (a black hole-torus system or hypermassive NS), 
which will lose mass through neutrino-driven and mag- 
netohydrodynamical outflows. The corresponding mat- 
ter will increase the ejecta mass and will ultimately have 
to be taken into account for reliable predictions of the 
properties of electromagnetic counterparts of NS merg- 
ers. The robustness of the nucleosynthesis outcome has 
to be explored concerning variations connected to un- 
certainties of the nuclear reaction rates. Further work 
is also needed to address how NS mergers as r-process 
sources fit into chemical evolution scenarios of the Milky 
Way, which should explain the observations of r-element 
enhanced metal-poor stars. Finally, the capabilities of 
various observational facilities have to be evaluated in 
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view of the bounds on the observable features set by our 
survey. 
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